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Abstract 

The electric dipole form factor (EDFF) of the nucleon stemming from the QCD 
term and from the quark color-electric dipole moments is calculated in chiral 
perturbation theory to sub-leading order. This is the lowest order in which the 
isoscalar EDFF receives a calculable, non-analytic contribution from the pion cloud. 
In the case of the 6 term, the expected lower bound on the deuteron electric dipole 
moment is \dd\ ^ 1.4-10 _4 # e fm. The momentum dependence of the isovector EDFF 
is proportional to a non-derivative time-reversal-violating pion-nucleon coupling, 
and the scale for momentum variation — appearing, in particular, in the radius of 
the form factor — is the pion mass. 



The electric dipole form factor (EDFF) completely specifies the parity (P) and time- 
reversal (T) -violating coupling of a spin 1/2 particle to a single photon [U [2] • At zero 
momentum, it reduces to the electric dipole moment (EDM), and its radius provides 
a contribution to the Schiff moment (SM) of a bound state containing the particle [3]. 
The full momentum dependence of the form factor can be used in lattice simulations to 
extract the EDM by extrapolation from a finite-momentum calculation pE] (in addition to 
the required extrapolations in quark masses and volume [5]). 

There has been some recent interest [JJ, [21 [61 El |9] in the nucleon EDFF motivated 
by prospects of experiments that aim to improve the current bound on the neutron EDM, 
\d n \ < 2.9 • lCT 13 e fm [lOj, by nearly two orders of magnitude [IT], and to constrain 
the proton and deuteron EDMs at similar levels [12]. We would like to understand the 
implications of a possible signal in these measurements to the sources of T violation at the 
quark level, which include, in order of increasing dimension, the QCD 9 term, the quark 
color-EDM (qCEDM) and EDM, the gluon color-EDM, etc. [131 EH]- Unfortunately, as 
with other low-energy observables, both the EDM and the SM of hadrons and nuclei are 
difficult to calculate directly in QCD. However, long-range contributions from pions can, 
to some extent, be calculated using the low-energy effective field theory of QCD, chiral 
perturbation theory (ChPT) JT5J [161 02] • ChPT affords a systematic expansion of low- 
energy observables in powers of Q/Mq C d, where Q represents low-energy scales such as 
external momenta and the pion mass m^, and Mqqd ~ 1 GeV denotes the characteristic 
QCD scale. (For introductions, see for example Refs. [TBI IT9] .) 

In Refs. [U [9] the nucleon EDFF stemming from T- violation sources of effective dimen- 
sion up to 6 was considered in ChPT to the lowest order where momentum dependence 
appears. It was argued [9] that the nucleon EDFF partially reflects the source of T vio- 
lation at the quark level. The various sources differ in particular in the expectation for 
the behavior of the isoscalar EDFF. For 9 and qCEDM, the isoscalar momentum depen- 
dence appears only at NLO. The nucleon EDFF from 9 was calculated at LO in Ref. [JJ, 
generalizing to finite momenta earlier calculations of the EDM [201 [21] and SM [3] . At 
this order, the momentum dependence is isovector and completely due to a T-violating 
coupling of the pion cloud to the nucleon, with a radius fixed by m\ [3]. In Refs. [21 [8] 
the EDFF calculation was extended to NLO, and corrections found to be significant. For 
the qCEDM, the nucleon EDFF was calculated at LO [9], at which order it is identical 
to that from 9. For the other sources of effective dimension 6, the quark EDM and the 
gluon color-EDM, the nucleon EDFF, including its isoscalar component, was calculated 
to NNLO [9] and found to be mostly determined by short- distance physics. 

Since the proposed deuteron experiment will probe the isoscalar combination of neu- 
tron and proton EDMs (in addition to T-violating two-nucleon effects), we present here 
results for the nucleon EDFF to NLO from both 9 and qCEDM, using SU(2) L x SU(2) R 
heavy-baryon ChPT [16]. For 9, we extend the calculation of Ref. [2] and reproduce 
the EDM results of Ref. [8], the latter obtained from a relativistic version of large- N c 
U(3)l x U(3)r ChPT, except for isospin-violating terms neglected in Ref. [8]. At this 
order, the isoscalar momentum dependence, and so the SM, is entirely due to an isospin- 
breaking term related to the nucleon mass splitting. As we are going to see, no new 
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undetermined parameters appear, other NLO contributions being given by non-analytic 
recoil corrections proportional to m n /m^r, where is the nucleon mass, and by another 
isospin-breaking term, related to the pion mass splitting. We use the non-analytic con- 
tributions to the isoscalar EDFF to provide an estimate of the minimum expected size 
of the deuteron EDM. The EDFF from the qCEDM depends at NLO on an additional 
T-violating pion-nucleon coupling, although it is unlikely that the difference could be 
isolated experimentally. 

For simplicity we focus here on QCD with two light quark flavors u and d, most relevant 
for low momenta Q ~ m^, and consider as explicit degrees of freedom only nucleons, 
pions, and photons. In the framework of ChPT, the most general effective Lagrangian is 
built up using QCD symmetries as a guide, in particular the chiral SUl{2) x SU(2)r ~ 
SO (4) symmetry, which is spontaneously broken down to SU(2)l+r ~ 5*0(3). A power- 
counting argument must be used to order interactions according to the expected size 
of their contributions. In order to fulfill chiral-symmetry requirements, pions couple 
derivatively in the chiral limit, which brings to amplitudes powers of pion momenta. 
Chiral-symmetry-breaking terms involve the quark masses m u and m^, so they bring 
into amplitudes powers of the pion mass. Since nucleons are non-relativistic, we remove 
the large, inert nucleon mass from nucleon fields [16]. This gives one a chiral index 
(A) with which to order terms in the Lagrangian [ISlinj, i.e. £= Ea £(A) - For stron g 
interactions, the index is given by A = d + n/2 — 2, where n is the number of fermion fields 
and d counts the number of derivatives and powers of the pion mass. Electromagnetic 
interactions are proportional to the small proton charge e = y/Aira em , and it is convenient 
to account for factors of e by enlarging the definition of d accordingly. Each interaction 
is associated with a parameter, or low-energy constant (LEC), which can be estimated 
using naive dimensional analysis (NDA) [22j [13]. In this case, the index A tracks the 
number of inverse powers of Mqqd ~ 2nF n ~ 1.2 GeV, with F n ~ 186 MeV the pion 
decay constant, associated with an interaction. (Note that since NDA associates the 
LECs of chiral-invariant operators with g s /47r, for consistency one should take the strong- 
interaction coupling g s ~ Air.) 

The theory can be enlarged in a straightforward way to include the delta isobar. Note 
that the delta isobar does not contribute to the nucleon EDFF at the order in which we 
work. As is the case for the nucleon anapole form factor [23], the structure of the delta 
interactions that would contribute at NLO vanishes in ChPT. The first nonvanishing delta 
contribution occurs at a higher order than we are considering here. 

The T-conserving terms that we will need consist of the following [191 12H 125], [2]: 

4V = ■ - ^tt 2 + Niv ■ VN - 2 jr(D^) ■ NrS^N (1) 
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Here 7r denotes the pion field in a stereographic projection of SO(4)/ SO(3), with D = 
1 + 7r 2 /F 2 ; N = (p n) T is a heavy-nucleon field of velocity and spin (S 71 = (0, c?/2) 
in the nucleon rest frame where t> M = (1,0)); and is the photon field. In addition, 
(Dfj,)ab = D~ 1 (5abd t j, — ee^abA^) is the pion covariant derivative; = <9 M + ir ■ (tt x 
D^iz)/F% — ieA^l + r 3 )/2 is the nucleon covariant derivative; and F^ u = d^A u — d u A /Ji is 
the electromagnetic field strength. The component of V 1 perpendicular to v M is written 

V*t = V - v^v ■ V, (3) 

and we use the shorthand notation 

= V ± £> f/ \ = r i T) t " ± D+^ri, (4) 

where T>* is defined through NT)* = XW. 

The pion-nucleon coupling g^ in Eq. ([1]) and the anomalous magnetic photon-nucleon 
couplings kq and n\ in Eq. (j2J) are not determined by symmetry, but are expected to 
be 0(1), and indeed g& = 1.3, = —0.12, and /ci = 3.7. The pion mass term in 
Eq. ([I]) originates in explicit chiral-symmetry breaking by the average quark mass fh = 
(ma + m u )/2; from NDA, m 2 = 0{w,Mqcd)- The contribution to the nucleon mass from 
a similar term, the nucleon sigma term, has been removed by an appropriate definition 
of the heavy-nucleon field, and the surviving nucleon-pion interactions do not contribute 
below. The Goldberger-Treiman relation g^NN = 2gxm N / 'F n holds in the two lowest 
orders, and a term in provides an 0(m%/MQ CD ) correction that accounts [25] for 

the so-called Goldberger-Treiman discrepancy. 

In Eq. ([2]) we include explicitly the leading isospin-breaking interactions [25] stemming 
from the quark mass difference, rrid — m u = 2em, and from short-range electromagnetic 
effects. The pion mass splitting is dominated by the electromagnetic contribution 5m 2 = 
0(a ern MQ CD /4:7i)] because this is, numerically, of 0{em\j Mqcd) we book this term in 

^str/em- The quark-mass contribution to the pion splitting is suppressed by a further 
em n / 'Mqcd, one order down in the expansion. Thus, to the accuracy in which we work 
here, <5m 2 is the observed pion mass splitting, 5m 2 = 1260 MeV 2 |26j. Note that with the 
way we have written the splitting, in this paper m v stands for the neutral pion mass, the 
charged pion mass squared being m 2 ± = m 2 + bm\. 

The quark-mass contribution to the nucleon mass difference, 5m^, is expected to be 
0(eml/M Q c D ) and it is evaluated to be 5m N = 2.26 ± 0.57 ± 0.42 ± 0.10 MeV from 
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lattice simulations [27], which is in agreement with an extraction from charge-symmetry 
breaking in the pn — > dn° reaction [28]. This is consistent with the NDA expectation 
that the corresponding electromagnetic contribution, 5m N , is 0{a em MQ C D/^) and thus 
somewhat smaller; using the Cottingham sum rule [29], SrriN = — (0.76 ± 0.30) MeV. In 

r (2) 

our power counting, <5mjv appears only in C-J r j em - It gives rise to different multi-pion 
interactions than does the Sm^ term, but these multi-pion interactions do not matter 
below. Thus, if desired, the dominant effect from Sm^ can be incorporated in our results 
with the simple substitution 5m N —> 5m N + 5m N . The isosp in-breaking nucleon mass 
term can be eliminated with simultaneous redefinitions of the pion and nucleon fields [30] . 
at the expense of new interactions. In the present case, the relevant terms amount to 
adding to Eq. (T2]) 

^% em = -5m N (nxv.Dn) 3 - S -^NT 3 N. (5) 

As a check, we have performed the calculation with and without Eq. ([5]), finding the same 
result, as required by field-redefinition invariance. 

We consider two sources of T violation at the scale Mqcd'- the QCD 9 term and the 
qCEDM. In terms of a gluon field strength G°L, and an appropriate choice [3TJ[32] of quark 
fields q = {u d) T , we can write them as 

Cf D = mj qi lb q - % -q (d + J 3 r 3 ) cr^ l5 X a q G% (6) 

where 

m, = ^^ = ^(l-e 2 ) = 0(m) (7) 
m u + m d 2 

and do (ofa) is the isoscalar (isovector) qCEDM. The first term in Eq. (j6]) represents the 
effect of the 9 term under the assumption that 9 is small, as inferred from the bound [H)] 
on the neutron EDM. (For the more general case, see Ref. [32].) The second term in 
Eq. ([6]) is the QCD manifestation of sources of T violation at a high scale Mp. At the 
Standard Model scale it is represented by dimension-6 operators containing the media- 
tor of electroweak-symmetry breaking [131 EH] , which at lower scales picks up a vacuum 
expectation value that can be traded for fa. We write [H] 

in terms of a dimensionless factor 5. The size of 5 depends on the exact mechanisms of 
electroweak and T breaking and on the running to the low energies where non-perturbative 
QCD effects take over. The minimal assumption is that it is 0(g s /4:n), with g s the strong- 
coupling constant, but it can be much smaller (when parameters encoding T- violation 
beyond the Standard Model are small) or much larger (since the first-generation Yukawa 
couplings are unnaturally small). 
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The implications of T violation to low-energy observables depend on the way Eq. OS) 
breaks other QCD symmetries, in particular chiral symmetry [2J [321 [S3]- The 9 term is 
the fourth component of the same SO (A) vector P = (qrq, qi^^q) that leads to isospin 
breaking [TJ [2j [32] . Therefore, it generates EFT interactions that transform as T- violating, 
fourth components of SO (4) vectors made out of hadronic fields, with coefficients related 
to those of T-conserving interactions. Similarly, the qCEDM breaks chiral symmetry 
as a combination of fourth and third components of two other SO (4) vectors [21 [9j 133] . 
As in the T-conserving case, we can use NDA to estimate the strength of the effective 
interactions, and continue to label terms in the effective chiral Lagrangian by the powers 
of Mqq D . Details of the construction of the Lagrangian from these terms are discussed 
in Refs. [321 [33]. 

The relevant terms here are the pion-nucleon interactions 



(<<) 



N (g ir ■ t + 5-i7r 3 ) N 



(9) 



and 



( n+ l) z 

L T - fib 



(D^) ■ N (h 7z + M 3 r) S»N 



F W D I" 10 FID 
and the short-range contributions to the nucleon EDM 



(10) 



r (n+2) r (n+3) 



2* { (l - Jp) + 4 (r 3 - ^ (TrVa - 7r 3 7r ■ r)) 



FID' 



Here hi, di, and d[ are parameters of sizes 



N 



go 



o 



m; 



Mi 



~mlM Q cD 



QCD 



5i 



hn = 9 



1 1 




and 



(12) 



(13) 



d ,i,d[ = O 



m„ 



M qcd 



(14) 



The isoscalar (do) and isovector (d%, d[) contributions to the nucleon EDM occur for both 
T-violation sources. Direct short-range contributions to the momentum dependence of 



the EDFF first appear in , being further suppressed by 0(Q / Mq C d)- For 9, 
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and among the T-violating irN interactions only the 1 = interactions with coefficients 
g and ho appear [JJ. Because of the link with isospin violation 



mJm N - 5m N - 
9o = z — d^——d. 15 
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A similar connection exists between ho and the leading isospin breaking in the pion- 
nucleon vertex. For qCEDM, n = —1 and all ttN terms should be included. Since in this 
case there is no analogous link to T-conserving quantities, one cannot do better than the 
NDA estimates f[T2j) and (TX5T) without lattice or dynamical- model input. 
The T-violating current-current nucleon-electron interaction is of the form 

iT = -iee(r)Ye(l) D^ u (q) N(p')r ed (q,K)N(p), (16) 

where e(l) (N(p)) is an electron (nucleon) spinor with momentum I (p) and D^q) = 
—iiilnu/q 2 + • • •) is the photon propagator with q 2 = (p — p') 2 = —Q 2 < 0. The nucleon 
electric dipole current J^ d can be expressed in terms of q = p — p' and K = (p + p')/2 as 
an expansion in powers of Qjm^ that reads [JJ [32j [9] 



J? d (q,K) = 2{F (Q 2 ) + F 1 (Q 2 )T 3 ) 
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where Fq(Q 2 ) (Fi(Q 2 )) is the isoscalar (isovector) EDFF of the nucleon. We will write 

F i (Q 2 ) = d i -S' i Q 2 + H i (Q 2 ), (18) 

where di is the EDM, the SM, and Hi(Q 2 ) accounts for the remaining Q 2 dependence. 

The form factor itself can be expanded in powers of Q/Mqcd- The leading-order (LO) 
contributions to the current, which are 0(egoQ/ (2nF n ) 2 ), have been calculated in Refs. 
[U [9]. They include the unknown short-range contributions in £^™ +2 ' ) ffTTj) . and loop 
diagrams made out of T-violating interactions in OH]) and T-conserving interactions 
in £g°L em (DO)- Here we focus on the next-to-leading order (NLO), that is, terms of relative 
0(Q/Mqcd)- They are made of diagrams with one insertion of interactions in £^? +3 ' ) ( TTTj) . 
£^ ra+1 ' ) f lTUj) . or £^t r / em (E]). There are no new, unknown short-range parameters appearing 
at tree level, the recoil corrections in £^? +3 ^ ( TTTj) simply ensuring — together with those in 

^str/em © — the form ( FlTl) of the current. The loop diagrams contributing to the nucleon 
EDFF in NLO are shown in Figs. [TJ and [2J classified according to the combination of 
couplings that they contain. All other contributions to the EDFF are formally of higher 
order: they come from more powers of momenta in diagrams with the same number of 
loops, or from extra loops. 

The NLO diagrams of Fig. [Hare built from the leading interactions in Eqs. ([1]) and (J9j), 
plus one insertion of an operator from Eq. (|2]). Diagrams [lja,b,c) represent a correction 
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Figure 1: One-loop diagrams contributing to the nucleon electric dipole form factor in 
sub-leading order coming from one insertion of an £^/em °P era tor. Solid, dashed and 
wavy lines represent nucleons, pions and (virtual) photons, respectively; single filled circles 
stand for interactions from C } , while double circles for interactions from C } , and/or 

str/em str/em ' 

A£^y em ; squares represent the T- violating vertices from . For simplicity only one 
possible ordering is shown here. 



to the external energies, 



m N 

v-K = -J-(W^ T — > ( 20 ) 
2m N V 4 J ^ 2 ' v ' 

of a proton (— sign) or neutron (+ sign) in LO diagrams. (In the remaining NLO dia- 
grams, we set the right-hand side of these equations to zero.) Analogous insertions in the 
nucleon propagator are represented by diagrams [T](d,e,f). Diagrams [T](g,h,i) originate in 
the recoil correction in pion emission/absorption, while diagram QTj) arises from the mag- 
netic photon-nucleon interaction. Diagrams [TJ^k,l,m) represent an insertion of the pion 
(and nucleon, in the case of Eq. mass splitting in pion propagation, while diagram 
[U(n) contains the isospin- violating photon-pion interaction induced by the first term of Eq. 
(jSJ). These one-loop diagrams contribute to the current at order 0(egiQ 2 / (2iiF n ) 2 mN) . 
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Figure 2: Diagrams contributing to the nucleon electric dipole form factor in sub-leading 
order coming from one insertion of the T- violating vertex from C^ +l \ represented by a 
circled square. Other symbols are as in Fig. [TJ For simplicity only one possible ordering 
is shown here. 

The NLO diagrams in Fig. [5] are built from the leading interactions in Eq. ([T]) with 
one insertion of an operator from Eq. ( 1T0"|) . Diagrams 2(a,b) stem from the sub- leading 
pion-nucleon couplings /z 0i i, and diagrams 2(c,d,e,f) from the sub-leading coupling h 2 , 
present only for qCEDM. These one-loop diagrams contribute to the current at order 
0(ehiQ 2 /(2irF n ) 2 ), which is precisely the same order as the diagrams in Fig. [TJ 

The diagrams in Figs. [TJ and [2] can be evaluated in a straightforward way. We use 
regularization in d spacetime dimensions, and define 

L=^- 7s + ln47r, (21) 

where je = 0.557 ... is the Euler constant. The LO loop contributions depend on a renor- 
malization scale fi but this dependence is compensated for by the contribution from the 
short-range interactions in Eq. (fTTJ) . The NLO diagrams are finite in this regularization 
scheme. 

Most of the diagrams actually vanish when the on-shell conditions ( 1191) and ( 120]) are 
consistently enforced. Diagrams (a,b) in Fig. [2] vanish due to isospin. Since diagrams 
[2](c,d,e,f) vanish too, the EDFF to this order depends only on the leading T-violating 
parameters g ,i through Fig. [TJ Diagram [TJ(j) vanishes due to its spin structure and there- 
fore the EDFF does not depend on the anomalous magnetic moments, either. Diagram 
[TJ(h) gives both isoscalar and isovector contributions. The remaining non-vanishing dia- 
grams are[]Ja,d,k). Neglecting T-conserving isospin violation, these diagrams give purely 
isovector results. In the case of 9, the results are proportional to egAgo/(2TrF n ) 2 , as in 
LO [TJ, times the recoil suppression factor m^/mjq. For qCEDM, there is an additional 
momentum-independent contribution proportional to g~\. 
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We have checked each of the isospin-breaking contributions in two ways. The contri- 
butions from 5m% come through diagrams [TJk,l,m). Because the LO EDFF originates 
entirely in charged-pion diagrams, these contributions can be obtained alternatively by 
evaluating the LO EDFF with m^ + Sm^, then expanding in powers of 5m\jm\: 
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The resulting EDFF is thus isovector. Including the nucleon mass difference 5m^, di- 
agrams HJa,d) generate an additional isoscalar contribution. If we instead use the field 
redefinition of Ref. |30j, we drop the 5m^ term in Eq. ( 1201) and in internal nucleon lines. 
The same result is obtained from extra contributions oc 5m^ in diagrams [T^k,l,m), and a 
new isospin-breaking photon-pion coupling, depicted in diagram QJn). 

The diagrams in Fig. [1] contribute to both isoscalar and isovector EDMs. Taking the 
NLO contributions together with the LO from Refs. [TJ [5], we have 
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The LO piece in Eq. f [24l) . which depends on go and is non-analytic in m^, is, with 
the use of the Goldberger-Treiman relation, the result of Ref. [20], which holds also for 
the qCEDM [9J. The short-range isovector combination d\ + d[ absorbs the divergence 
and fj, dependence of the LO loop. The short- and long-range contributions to the EDM 
are in general of the same size, but a cancellation is unlikely due to the non-analytic 
dependence on m n of the pion contribution. The isoscalar parameter do is not needed for 
renormalization at this order, but there is no apparent reason to assume its size to be 
much smaller than NDA either. 

At NLO, the EDM receives finite non-analytic corrections, which depend also on cji 
for qCEDM. From Eqs. (F23D an d (ED we see that, as usual in baryon ChPT, the NLO 
contributions are enhanced by it over NDA. However, the other dimensionless factors are 
not large enough to overcome the m- K jm N suppression. Setting \i to m N as a representative 
value for the size of di [20], the NLO term in Eq. (Eq. flSJ) is about 15% (10%) of the 
leading non-analytic term in Eq. (|24|) . indicating good convergence of the chiral expansion. 
The isovector character of the LO non-analytic terms is approximately preserved at NLO. 
Isospin-breaking contributions, although formally NLO, are pretty small, amounting to 
15-20% of the total NLO contribution. 

In the case of 9 we can use Eq. f fl5|) and expect 
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for the neutron EDM and 
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for the proton EDM, using the lattice QCD value 5m N /2e = 2.8 MeV [27]. Non-analytic 
NLO corrections are therefore somewhat larger for the proton EDM, but this difference 
is unlikely to be significant in light of our ignorance about the size of short-range contri- 
butions. 

The non-analytic terms in Eq. fT23|) represent a lower-bound estimate for the size of the 
nucleon isoscalar EDM, as the short-range contribution d is nominally of lower order. The 
expected lower bound on the nucleon isoscalar EDM might have implications for proposed 
experiments on EDMs of light nuclei. In these cases, there will be additional many-nucleon 
contributions, but the average of the one-nucleon contributions still provides an estimate 
of the order of magnitude of the expected nuclear EDM. For the deuteron, the average 
one-nucleon contribution is exactly do and, in the case of 9, we expect for the deuteron 
EDM 
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(27) 



Therefore, if there are no cancellations, a deuteron EDM signal from 9 is expected to be 
larger than about 10% of the neutron EDM signal. 

Note that short- and long-range physics cannot be separated with a measurement of the 
neutron and proton EDMs alone. On the other hand, the momentum dependence of the 
EDFF is completely determined, to the order we are working, by long-range contributions 
generated by go- It is therefore the same for 9 and qCEDM. It turns out that the isoscalar 
form factor receives momentum dependence only from isospin-breaking terms, while there 
is a non-vanishing correction to the isovector momentum dependence also from isospin- 
conserving terms. 

The variation of the form factor with Q 2 can be characterized at very small momenta 
by the electromagnetic contribution to the nucleon SM, the leading and sub-leading con- 
tributions of which we find to be 



S' 



egAQo TT 5m 



N 



Q(2nF 7r ) 2 ml2 m w 



egAQo 



6(2nF 7T ) 2 ml 



57r m^ 5m 2 
4 m N 



mi 



(28) 
(29) 



The LO, isovector term is the result of Refs. (3j [T]. While the EDM vanishes in the 
chiral limit, the isovector SM is finite. The NLO correction, which agrees with the 9 
result of Ref. [8] when T-conserving isospin violation is neglected, vanishes in the chiral 
limit but gives a relatively large correction to the isovector SM of about 60%, due to the 
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numerical factor 57r/4. Again, the isospin-breaking corrections are relatively small, and, 
as a consequence, at NLO the SM remains mostly isovector. 

To this order, the SM is entirely given, apart for g , by quantities that can be deter- 
mined from other processes. In the case of 6, we can again use Eq. (115)) to estimate 



S' 



egA 



ir(5m N ) 2 - 



12(2ttF^) 2 2em\ 



-5.0 • l(n 6 # efm d 



egA 



5m 



N 



12(2ttF^) 2 em\ 



1 - 



5tt m n 
4 rriN 



5ml 



mi 



9 ~ 6.8 • 10" 5 # efm J 



(30) 
(31) 



where again we used the lattice-QCD value [27] for 5m^j1e. From these results we can 
straightforwardly obtain the SM for the proton and the neutron. Although we could again 
use the isoscalar component as an estimate for a lower bound on the deuteron SM, there 
could be potentially significant contributions from the deuteron binding momentum. 

The full momentum dependence of the EDFF is given in addition by the functions 
Hi(Q 2 ) introduced in Eq. (jig]) . 



H (Q 2 ) 
Hi(Q 2 ) 



4eg A g 3tc 5m N (1 ) / Q 2 \ 



15(2vrF^) 2 4 m n 



4eg A go 



15(2tt^) 2 
Here, the LO term, 

,(o) 



(0) 



4m 2 / 



4m 2 / ' 
7vr m^ h{1) 

8 TTljV 1 



Ami 



2Sm l m ( Q 



mi 



Am 2 



(32) 



(33) 



hr(x) 



15 

T 



In 



y/l + l/x + 1 



2 1 



y/l + 1/X- 1 

is the one calculated in Refs. [1] |9], while we now obtain the NLO isovector functions 



(34) 



h^(x) = ~ 3(1 + 2x)h ( u'{x) - lOx 



Mi 



7 



and 



h?(x) 

and the NLO isoscalar function 



4(1 + x) 



hf\x)-hx< 



hjp (x) = 5 ( —= arctan y/x — 1 + — j . 



(35) 
(36) 

(37) 
x 2 + 0(x 3 ) 



In compliance with the definition of if,, the four functions behave as h\ '{x) 
for x <C 1. 

As in lowest order, the momentum dependence is fixed by the pion cloud. Thus the 
scale for momentum variation is determined by 2m n . Both the SM and the functions 
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#0,1 



15 (2nF„) 2 
2.0 \ 



/ 



0.5 



1.5 



1.0 




Q 1 



0.5 



1.0 



13 



-2.0 Ami 



Figure 3: The isovector Hi(Q 2 ) in LO (dashed line) and LO+NLO (solid line), and the 
isoscalar H (Q 2 ) in LO+NLO (dash-dotted line), both in units of 4egA#o/15 (Z^F^) 2 , as 
functions of Q 2 , in units of 4m 2 . 

H 0il (Q 2 ) are testable predictions of ChPT. In Fig. |3]we plot the LO hf ] , the LO+NLO 

combination — (T^m^/Sm^M — (25m^./m^)^ 1 ' ) , and the NLO — (3ir5mN / Am n )h^ 
as functions of Q 2 . We use the same values of parameters as before. As for the SM, NLO 
corrections can be significant, but the isosp in-breaking contributions are small. 

In summary, we have calculated the nucleon electric dipole form factor due to the 9 
term and to the quark color-electric dipole moment in sub-leading order in ChPT, includ- 
ing isospin-breaking effects. The chiral expansion seems to be converging, although NLO 
corrections are enhanced by extra factors of tt. Under the assumption that higher-order 
results are not afflicted by anomalously-large dimensionless factors, the relative error of 
our results at momentum Q should be ~ (Q/Mqcd) 2 - The NLO isospin-breaking con- 
tributions are relatively small and could be overcome by isospin-conserving contributions 
at NNLO. We have shown that at NLO the EDFF includes both isoscalar and isovector 
components, with a Q 2 dependence determined by non-derivative T- violating pion-nucleon 
couplings and the pion mass. The isoscalar momentum dependence is entirely due to the 
nucleon mass splitting. We have provided a lower-bound estimate for the isoscalar nucleon 
EDM, expected to set also the minimum size of the deuteron EDM. A full calculation of 
the latter in ChPT can now be performed. 
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